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Abstract 


This paper is a proceedings version of [6], in which we state a Quantum Ergodicity (QE) 
theorem on a 3D contact manifold, and in which we establish some properties of the Quantum 
Limits (QL). 

We consider a sub-Riemannian (sR) metric on a compact 3D manifold with an oriented 
contact distribution. There exists a privileged choice of the contact form, with an associated 
Reeb vector field and a canonical volume form that coincides with the Popp measure. We state 
a QE theorem for the eigenfunctions of any associated sR Laplacian, under the assumption 
that the Reeb flow is ergodic. The limit measure is given by the normalized canonical contact 
measure. To our knowledge, this is the first extension of the usual Schnirelman theorem to 
a hypoelliptic operator. We provide as well a decomposition result of QL’s, which is valid 
without any ergodicity assumption. We explain the main steps of the proof, and we discuss 
possible extensions to other sR geometries. 

1 Introduction 

The property of Quantum Ergodicity (QE) is defined as follows. Let M be a metric space, endowed 
with a measure fi defined on a suitable compactification of M, and let T be a self-adjoint operator 
on L^(M,/i), bounded below and having a compact resolvent (and hence a discrete spectrum). 
Let ((/>ra)ngN* be a (complex-valued) Hilbert basis of L^{M,fi), consisting of eigenfunctions of T, 
associated with the ordered sequence of eigenvalues Ai ^ ^ A„ ^ . We say that QE holds 

for T is there exist a probability measure v on M and a density-one sequence (nj)jgN* of positive 
integers such that the sequence of probability measures converges weakly to v. The 

measure v may be different from some scalar multiple of /r, and may even be singular with respect 
to /i. 

Microlocal versions of QE hold true in general and are stated in terms of pseudo-differential 
operators, yielding a result of the kind 
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for every classical symbol a of order 0, where 1/ is a lift to T*M of the measure u and S is a subset 
of T*M (that have to be identified). 

Such a QE property provides an insight on the way eigenfunctions concentrate when considering 
highfrequencies. 

In the existing literature, QE results exist in different contexts, but always for elliptic operators. 
The first historical QE result is due to Shnirelman in 1974 (see [20] for a sketch of proof), stating 
that, on a compact Riemannian manifold {M,g), if the geodesic flow is ergodic, then we have QE 
for any orthonormal basis of eigenfunctions of the Laplace-Beltrami operator, with dv equal to the 
normalized Riemannian volume (and diy equal to the normalized Liouville measure). The complete 
proof has been obtained later in [5, 21]. Then, the Shnirelman theorem has been extended to the 
case of manifolds with boundary in [9, 23], to the case of discontinuous metrics in [13] and to the 
semiclassical regime in [10]. 

We provide here a Shnirelman theorem in a hypoelliptic context. 

Theorem 1. Let M be compact connected 3D manifold (without boundary), endowed with an 
arbitrary volume form dp, and with a Riemannian contact structure. We consider the associated 
sub-Riemannian Laplacian As_r. If the Reeb flow is ergodic for the Popp measure, then we have 
QE, and the limit measure is the probability Popp measure. 

Several definitions are due to understand the content of this theorem. In Section 2, we are thus 
going to recall: 

• what is a sub-Riemannian Laplacian; 

• what are a contact structure, the Reeb flow, and the Popp measure. 

We will then provide a more general (microlocal) version of Theorem 1. We also provide a second 
result, valid without any ergodicity assumption, stating that every quantum limit that is supported 
on the characteristic manifold of the sR Laplacian is invariant under the Reeb flow, and that most 
QL satisfy the previous assumption on their support. 

We will explain the main steps of the proof in Section 3, and we will conclude in Section 4 by 
giving some perspectives and open problems. 


2 Sub-Riemannian Laplacians in the 3D contact case 

2.1 General definition 

We first recall that a sub-Riemannian (sR) structure is a triple {M,D,g), where M is a smooth 
manifold, U is a smooth subbundle of TM (called horizontal distribution), and g is a fibered 
Riemannian metric on D. 

Example 1. We speak of a contact structure when M is of dimension 3 and D is a contact 
distribution, that is, we can write D = kera locally around any point, for some one-form a such 
that a f\ da fQ (locally). The local Darboux model in at the origin is given hy a = dz x dy, 
and then, locally, we can write D = Span(X, Y), with X and Y the vector fields defined hy X = dx 
and Y = dy — xdz- Note that D is of codimension one, and the Lie bracket [X, Y] = —dz generates 
the missing direction. 
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Sub-Riemannian Laplacian. In order to define a sub-Riemannian Laplacian Asr, let us choose 
a smooth volume form dfi on M, the associated measure being denoted by /i. Let fi) be the 

set of complex-valued functions u such that |up is ^-integrable over M. Then, in whole generality, 
—Asr is the nonnegative self-adjoint operator on defined as the Friedrichs extension of 

the Dirichlet integral 

QW = f IM<^llg- dfi, 

Jm 

where the norm of d(f> is calculated with respect to the co-metric g* on T*M associated with g. 
The sR Laplacian Asr depends on the choice of g and of dg. 

Let us give two equilavent definitions of Asr: 

• We have 

AsRcj) = div^(Vsi?0) 

for every smooth function (/) on M, where div^ is the divergence operator associated with 
the volume form dg, defined by Cxdg = div^(X) dg for any vector field X on M, and where 
the horizontal gradient Xsr4> of a smooth function cj) is the unique section of D such that 
gqi'^sR(l>iq),v) = d(j){q).v, for every v G Dq. 

• If (Xi,..., Xm) is a local t^-orthonormal frame of D, then Vsr4> = and 


= - ^ X*X, = Y, + cliV;,(X,)X,) 

i=l i=l 


Hypoellipticity. A well known theorem due to Hormander (see [11]) states that, under the 
assumption Lie(Z?) = TM, the operator Asr is hypoelliptic (and even subelliptic), has a compact 
resolvent and thus a discrete spectrum Ai ^ A 2 ^ ^ A„ —>■ -foo. 

Remark 1. Denoting by hx{q,p) = {p,X(q)) the Hamiltonian associated with a vector field X 
(in canonical coordinates), the principal symbol of —Asr is 

m 

ap{—AsR) = = g*, 

i=l 

(it coincides with the co-metric g*), and its sub-principal symbol is zero. 

Actually, in our main result, Asr may be any self-adjoint second-order differential operator 
whose principal symbol is g*, whose sub-principal symbol vanishes. 

Remark 2. The Hamiltonian function Ag* induces a Hamiltonian vector field Xg which generates 
the so-called normal geodesics of the sR structure (see [17]). 

2.2 The 3D contact case 

We assume that the manifold M is compact and of dimension 3, and that D is an oriented contact 
distribution, that is, we can write D = ker a for some one-form a globally defined such that 
a f\ da (see [ 6 ] for the extension to the non-orientable case). 

In order to define the canonical contact measure and the Reeb vector field, we need to normalize 
the contact form defining D. There exists a unique contact form ag such that dag{X,Y) = 1 for 
any positive g-orthonormal local frame {X,Y) of D. Equivalently, (dag)|^ coincides with the 
volume form induced by 5 on D. 
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Popp measure. We define the density dP = \ag A dag\ on M. In general, dP differs from 
d^. The corresponding measure P is called the Popp measure in the existing literature (see [17], 
where it is defined in the general equiregular case). Of course, here, this measure is the canonical 
contact measure associated with the normalized contact form ag. We also define the probability 
Popp measure 

P 

^ ~ P{M)' 

Reeb flow. The Reeb vector field Z of the contact form ag is defined as the unique vector field 
such that Lzctg = 1 and iz dag = 0. Equivalently, Z is the unique vector field such that 

[X, Z] e D, [Y, Z] e £), [X, y] = -Z mod D, 

for any positive orthonormal local frame {X^Y) of D. The Reeb flow TZt is defined as the flow 
generated on M by the vector field Z. 

Remark 3. The measure u is invariant under the Reeb flow TZf This invariance property is 
important in view of identifying a QE result. 

Remark 4. Denoting by hz the Hamiltonian associated with Z, geodesics with highfrequencies 
in hz oscillate around the trajectories of ±Z. Erom the point of view of semiclassical analysis, this 
part of the dynamics is expected to be the dominant one. 

2.3 The main results 

We consider the cotangent space (T*M,w), endowed with its canonical symplectic form. Let 
E C T*M be the characteristic manifold of —/\sR- We have E = (5*)~^(0) = (annihilator of 
D), and E coincides with the cone of contact forms defining D. In particular, we have 

E = {(g, sag{q)) G T*M \ q G M,s G M}. 

In the 3D contact case, we have as well E = {hx = hy = 0} \i D = Span(X,F), with {X,Y) a 
local ( 7 -orthonormal frame of D. 

Shnirelman theorem in the 3D contact case. We are now in a position to provide a more 
precise statement of Theorem 1. 

Theorem 2. Let M be a smooth connected compact three-dimensional manifold, equipped with an 
arbitrary smooth volume form dp. Let D G TM be a smooth oriented contact subbundle, and let 
g be a smooth Riemannian metric on D. Let be the sR Laplacian associated with the contact 
sub-Riemannian structure {M,D,g) and with the volume form dp. 

We assume that the Reeb flow is ergodic^. Then, for any real-valued orthonormal Hilbert basis 
((/>ri)rtGN* of , p) consisting of eigenfunctions of Asr associated with the eigenvalues 
there exists a density-one sequence of positive integers such that 

lim ,(/)„ ) = i / {a{q,ag{q)) Y a{q,-ag{q))) dv, 

J^+oo Z 

for every classical pseudo-differential operator A of order 0 with principal symbol a. 

In particular, we have 

lim [ f\(j).„_fl^ dp = [ fdu, 
j^+oo JM Jm 

for every continuous function f on M. 

^The Reeb flow is ergodic on (M, u) if any measurable invariant subset of M is of measure 0 or 1. 
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Here, the notation ( , ) stands for the (Hermitian) inner product in 

We prove in [6] that the above result still holds with a complex-valued eigenbasis, provided the 
principal symbol a satisfies a{q,ag{q)) = a{q, —ag{q)), for every q € M. We also prove that the 
result is valid as well if D is not orientable, without any additional assumption (note that, in that 
case, the contact form and thus the Reeb vector held Z are not dehned globally). 

Note that the subsequence of density one depends on the choice of the eigenbasis ((/>„)„gn* • We 
stress that our result is valid for any choice of a smooth volume form on M. 

Remark 5. The classical Shnirelman theorem is established in the Riemannian setting under the 
assumption that the Riemannian geodesic how is ergodic on {S*M, Xl)-, where the limit measure is 
the Liouville measure Al on the unit cotangent bundle S*M of M. In contrast, here the Liouville 
measure on the unit bundle g* = 1 has inhnite total mass (where g* is the co-metric on T*M 
associated with g), and hence the QE property cannot be formulated in terms of the geodesic how. 

Another interesting difference is that, in the Riemannian setting, QE says that most eigenfunc¬ 
tions equidistribute in the phase space, whereas here, in the 3D contact case, they concentrate on 
E = £)-*-, the contact cone that is the characteristic manifold of Asr- 



Riemannian case 

SR case 

Ergodicity assumption 

geodesic how on (S'*M, Liouville) 

Reeb how on (M, Popp) 

Quantum limit 

Liouville measure on S*M 

Popp measure on M 

Microlocal concentration 

on S*M 

on STi = E n {\hz\ = 1} 


Remark 6 (Examples of ergodic Reeb hows in dimension 3). Let us give two general constructions 
providing examples of ergodic Reeb hows on 3D contact manifolds: 

• Geodesic hows: Let X be a two-dimensional compact Riemannian surface, endowed with a 
Riemannian metric h, and let M = S*X be the unit cotangent bundle of X. The closed 
three-dimensional manifold M is then naturally endowed with the contact form a dehned as 
the restriction to M of the Liouville 1-form X = pdq. Let Z be the associated Reeb vector 
held. Identifying the tangent and cotangent bundles of X with the metric h, the set M is 
viewed as the unit tangent bundle of X. Then Z is identihed with the vector held on the unit 
tangent bundle of X generating the geodesic how on S*X. Therefore, with this identihcation, 
the Reeb how is the geodesic how on M. 

This geodesic how is ergodic for instance if the curvature of X is negative. 

• Hamiltonian hows: Let (IE, ui) be a symplectic manifold of dimension 4, and let M be a 
submanifold of W of dimension 3, such that there exists a vector held t; on a neighborhood 
of M in W, satisfying LyUJ = uj (Liouville vector held), and transverse to M. Then the one- 
form a = ivOJ is a global contact form on M, and we have da = lo. Note that, if uj = dX is 
exact, then the vector held v dehned by LyUj = A is Liouville (in local symplectic coordinates 
{q,p) on W, we have v = pdp). If the manifold M is moreover a level set of an Hamiltonian 
function h on W, then the Reeb how on M (associated with a) is a reparametrization of the 
Hamiltonian how restricted to M. 

If D = ker a is moreover endowed with a Riemannian metric g, then ag = ha for some 
smooth function h (never vanishing). Let us then choose the metric g such that h = 1. Then, 
the Reeb how is ergodic on (M, v) if and only if the Hamiltonian how is ergodic on (W, w^). 
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Quantum limits in the 3D contact case. Let (V’j)jeN* be an arbitrary orthonormal family of 
We set (a) = {Op{a)iljj for every j € N*, and for every classical symbol a of order 
0. The measure fij is asymptotically positive, and any closure point (weak limit) of is a 

probability measure on the sphere bundle S*M, called a quantum limit (QL), or a semi-classical 
measure, associated with the family 

Theorem 1 says that, under the ergodicity assumption of the Reeb flow, the probability Popp 
measure v, which is invariant under the Reeb flow, is the “main” quantum limit associated with 
any eigenbasis. 

The following result provides an insight on quantum limits of eigenfunctions in the 3D contact 
case in greater generality, without any ergodicity assumption. In order to state it, we identify 
S*M = [T*M \ {0})/(0, +oo) with the union of the unit cotangent bundle U*M = {g* = 1} and 
of the sphere bundle SH = (S \ {0})/(0,+oo) which is a two-fold covering of M. Each fiber is 
obtained by compactifying a cylinder with two points at infinity. Moreover, the Reeb flow can be 
lifted to 5S. 

Theorem 3. Let (^n)neN* be an orthonormal Hilbert basis of L^[M,g), consisting of eigenfunc¬ 
tions of Asr associated with the eigenvalues (A„)„gN*- 

1. Let j3 be a QL associated with the family {4>n)n&i* ■ Using the above identification S*M = 
U*M U ST,, the measure (3 can be identified to the sum /3 = /3o + Poo of two mutually singular 
measures such that: 

• Po is supported on U*M and is invariant under the sR geodesic flow associated with the 
sR metric g, 

• Poo is supported on ST and is invariant under the lift of the Reeb flow. 

2. There exists a density-one sequence (nj)jgN* of positive integers such that, if P is a QL 

associated with the orthonormal family then the support of P is contained in ST, 

i.e., Po =0 in the previous decomposition. 

The proof of that result follows from arguments used to prove Theorem 1 (see [6] for details). 

2.4 The model example: compact Heisenberg flat case 

The simplest example is given by an invariant metric on a compact quotient of the Heisenberg 
group. The spectral decomposition of the Heisenberg Laplacian is then explicit and can serve as a 
model to derive our main result. 

Let G be the three-dimensional Heisenberg group defined as G = with the product rule 

(x, y, z) * (x', y', z') = {x + x',y + y', z^z' - xy'). 

The contact form an = dz -\- xdy and the vector fields Xh = dx and Yh = dy — xdz are left- 
invariant on G. Defining the discrete co-compact subgroup T of G by T = {(x, j/, z) G G | x,y € 
z G 27rZ}, we then define the three-dimensional compact manifold manifold Mr = T\G, 
and we consider the horizontal distribution Dr = kera//, endowed with the metric gR such that 
{Xr, Yr) is a gj^-orthonormal frame of Dr. With this choice, we have ctRg = aR. 

The Reeb vector field is given by Zr = —[Xr,Yr] = d^. The Lebesgue volume dg = dxdydz 
coincides with the Popp volume dP, and we consider the corresponding sub-Riemannian Laplacian 
Ar = Xfj -I- (here, the vector fields have divergence zero). 

We refer to this sub-Riemannian case as the Heisenberg flat case {Mr, Ar). The cotangent 
space T* Mr is endowed with its canonical symplectic form. 
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It is proved in [4] that the spectrum of —Ah is given by 


{^e,m — (2£ + 1)|to| I m G Z \ {0}, i G N} U — 27r(j^ + k^) | j, fc £ Z}, 


where is of multiplicity |m|. From this, we infer two things. 

The first (already known, see, e.g., [15]) is that the spectral counting function has the asymp¬ 
totics 


7V(A) ~ ^ 


A2 


^ ,2 _ ^(Mh) ,2 


1=0 


(21 + 1)2 8 


= -:-A^ = 


32 


as A —> -foo. It is interesting to compare that result with the corresponding result for a Riemannian 
Laplacian on a three-dimensional closed Riemannian manifold, for which JV(A) ~ C\^/^. 

The second thing is the following. Microlocally, in the cone Cc = > c(p^ -I- p^)} for 

some c > 0 (with local coordinates p = (px,Py,Pz) in the cotangent space), the sub-Riemannian 
Laplacian can be written as a commuting product 


—Ah = Rh^h = HhRh , 


with 

• Rh = ^Z'^Zh (pseudo-differential operator of order 1), 

• VIh = Uh + Vh (harmonic oscillator), where Uh and Vh are the pseudo-differential operators 
of order 1/2 defined by 

- op» (^) , ^ op- (^) . 

where the notation Op^ stands for the Weyl quantization. Note that [Uh,Vh] = ±id 
(according to the sign of Hzh) and exp(2i7rn//) = id. 

In terms of symbols, since Xh and Yh have divergence zero, we have factorized a(—AH) = 

+ K, (Ml syntol) as H^„| + (;fe) ) ■ 

3 The main steps of the proof of Theorem 2 

3.1 General concepts 

We fix an arbitrary (orthonormal) eigenbasis (</>n)neN* of Asr. We recall that the spectral counting 
function is defined by N(\) = #{n | A„ ^ A}. 

Definition 1. For every bounded linear operator A on L2(M,p), we define the Cesdro mean 
E(A) GCby 

^‘'‘> = xi;T=sjv?A) 

whenever this limit exists. We define the variance V (s4) G R"*" by 

I/(s4) = limsup-^ \{A(f>^,(f>ri)\^ . 

A^+oo tv (A) ^ 
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It is easy to see that V{A) ^ E{A*A), and that, if A is a compact operator on ^), then 

E{A) = 0 and V{A) = 0. 

Definition 2. The local Weyl measure is the probability measure on M defined by 

/ fdWA=Jim V / f\f>n\^dp, 

Jm a^-+oo tv (A) Jm 

for every continuous function / : M —>■ R, whenever the limit exists. In other words, 

n;A=weak liin V |(A„|V- 

A^+oo N{\) 


The microlocal Weyl measure VFa is the probability measure on S'^M = S{T*M), the co-sphere 
bundle, defined as follows: we identify positively homogeneous functions of degree 0 on T*M with 
functions on the bundle S*M. Then, for every symbol a : S*M —>■ R o/ order zero, we have 


a dW/\ 


IS*M 


1 


1 ™ 

A—>+00 W(A) 


A„<A 


where Op^ is a positive quantization, whenever the limit exists. In other words, Wa is the weak 
limit of the probability measures on S*M defined by 

^ (Op+(a)(?!>„, (?!>„). 


The microlocal Weyl measure, if it exists, does not depend on the choice of the quantization 
and of the choice of the orthonormal eigenbasis (because it is a trace). Moreover, Wa is even with 
respect to the canonical involution of S*M. 


General path towards QE. In order to establish QE, we follow the general path of proof 
described in [22]. Defining the spectral counting function by N{\) = ff{n \ A„ ^ A}, the first step 
consists in establishing a microlocal Weyl law: 

E{A)= [ adWA, (1) 

J ST* M 

for every classical pseudo-differential operator A of order 0 with a principal symbol a, where ST*M 
is the unit sphere bundle of T*M. This Cesaro convergence property can usually be established 
under weak assumptions, without any ergodicity property. 

The second step consists in proving the variance estimate 

V{A- E{A)iA)=Q. (2) 

The variance estimate usually follows by combining the microlocal Weyl law (1) with ergodicity 
properties of some associated classical dynamics and with an Egorov theorem. 

Then QE follows from the two properties above. Indeed, for a fixed pseudo-differential operator 
A of order 0, it follows from (2) and from a well known lemma^ due to Koopman and Von Neumann 
(see, e.g., [18, Chapter 2.6, Lemma 6.2]) that there exists a density-one sequence (n^ )jgN* of positive 
integers such that (Acfn^, (fuj) E{A) as j —>■ -l-oo. Using the fact that the space of symbols of 
order 0 admits a countable dense subset, QE is then established with a diagonal argument. 

^This lemma states that, given a bounded sequence (urtlneN of nonnegative real numbers, the Cesaro mean 
I converges to 0 if and only if there exists a subset S C N of density one such that (uAkeS converges to 

0. We recall that S is of density one if G S | fc < n — 1} converges to 1 as n tends to +oo. 







3.2 The microlocal Weyl law 

In order to prove Theorem 2, we follow the general path above, and we start by providing a 
microlocal Weyl law, identifying the microlocal Weyl measure in the 3D contact case. 

Theorem 4. Let A be a classical pseudo-differential operator of order 0 with principal symbol a. 
We have 


V +o(l)) [ {a{q,ag{q))-\-a{q,-ag{q)))dv, 

as A —>• +c». In particular, it follows that N{X) ~ and that 

E{A) = l-[ {a{q,ag{q)) + a[q,-ag{q))) dv = [ adW^, 

^ Jm iss 

where SH = (E \ {0})/(0, +oo) is the co-sphere bundle of S. 

Our proof of Theorem 4, done in [6], consists in first establishing a local Weyl law, by computing 
heat traces with the sR heat kernel and using the Karamata tauberian theorem, and then in infering 
from that local law, the microlocal one. The latter step is performed as follows. We first prove 
that, if the microlocal Weyl measure Wa exists, then Supp(Wa) C ST,. This follows from the fact 
that, outside of the characteristic manifold E, the operator is elliptic, and therefore classical 
arguments of [8, 12], using Fourier Integral Operators and wave propagation, yield that 

^ =0(A3/2), 

A„^A 

as A —^ + 00 , for every pseudo-differential operator of order 0 whose principal symbol vanishes in 
a neighborhood of E. Since we already know that iV(A) ~ CX^, the concentration on E follows. 
It can be noted that these arguments are as well valid for general sub-Riemannian structures, not 
only in the 3D contact case. 

Then, we prove that, since the horizontal distribution D is of codimension 1 in TM , and since 
we have already identified the local Weyl measure wa, then the microlocal Weyl measure Wa 
exists and is equal to half of the pullback of wa by the double covering ST —>■ M which is the 
restriction of the canonical projection of T*M onto M. This latter argument is also valid in a 
more general context, as soon as D is of codimension 1. 

Remark 7. A remark which is useful in order to understand why our result does not depend on the 
choice of the measure is the following. Let us consider two sR Laplacians and A^^ associated 
with two different volume forms (but with the same metric g). We assume that p 2 = h'^Li with 
h a positive smooth function on M. We define the isometric bijection J : if {M, 1x2) —t L^{AI,pi) 
by Jcf) = hcf). Then JA^^ = A^^ -f hAg_.^{h~^) id, and therefore A^^ is unitarily equivalent to 
A ^2 + W, where IF is a bounded operator. 

3.3 Birkhoff normal form 

Having in mind the model situation given by the Heisenberg flat case {Mh, Ah), in the general 
3D contact case we are able to establish a Birkhoff normal form, in the spirit of a result by Melrose 
in [14, Section 2], which implies in particular that, microlocally near the characteristic cone, all 
3D contact sub-Riemannian Laplacians (associated with different metrics and/or measures) are 
equivalent. 
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We define the positive conic submanifolds = {(g, sag{q)) € T*M \ ±s > 0} of T*M. In the 
Heisenberg fiat case, the characteristic cones are defined accordingly. 

Given /c G N U {+ 00 } and given a smooth function / on T*M, the notation / = Oi;(fc) means 
that / vanishes along E at order k (at least). The word flat is used when k = + 00 . 

Classical normal form. 

Theorem 5. Letq £ M be arbitrary. There exist a conic neighborhood Cq ofYi'^ in (T*M,u;) and a 
homogeneous symplectomorphism x from Cq to (T* MhjOjh), satisfying xiq) = 0 and x(E+nC'g) C 
E^, sueh that 

api-An) o X = (Tp{-Asr) + Oe(oo) 

Of course, a similar result can be given for E“. 

In order to establish this normal form, we first endow K® with a symplectic form uj, with 
an appropriate conic structure, and with an Hamiltonian function H 2 , such that, for any given 
contact structure and any q G M, there exists a homogeneous diffeomorphism xi from a conic 
neighborhood Cq of E+ to K® such that Xi^ = w + Oi;(l) and H 2 o xi = ap{—AsR). Thanks to 
the Darboux-Weinstein lemma, we modify xi info a homogeneous diffeomorphism X 2 such that 
X 2 W = w and H 2 o X 2 = o'p{—Asr) + Oi;(3). Finally, we improve the latter remainder to a flat 
remainder Oi;(oo), by solving an infinite number of cohomological equations in the symplectic 
conic manifold (]R®,ci)). 

Quantum normal form. By quantizing the above Birkhoff normal form, we obtain the following 
result. 

Theorem 6. For every q G M, there exists a (conic) microlocal neighborhood U ofTiq in T*M such 
that, considering all the following pseudo-differential operators as acting on functions microlocally 
supported in U, we have 

~Asr = Rfl + Vb + Oe(oo) 

where 

• Vb G Ik® is a self-adjoint pseudo-differential operator of order 0, 

• R and O are self-adjoint pseudo-differential operators of order 1, of respective principal sym¬ 
bols satisfying 


cpiR) — |^z| + Oi;(2), ap{n) — (h^x + ^f)/I^.z| + Oi;(3), 

• [R, n] = 0 mod , 

• exp(2i7rH) = id mod '1'“°°. 

The latter remainders are in the sense of pseudo-differential operators of order — 00 . In the flat 
Heisenberg case, there are no remainder terms, and we recover the operators Rr and Tip. The 
pseudo-differential operators R and H can be seen as appropriate perturbations of Rr and TIr, 
designed such that the last two items of Theorem 6 are satisfied. 

Note that the operators R and H depend on the microlocal neighborhood U under consideration. 
This neighborhood can then be understood as a chart in the manifold T*M, in which the quantum 
normal form is valid. 
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3.4 Sketch of proof of the variance estimate in the Heisenberg flat case 

As explained at the end of Section 3.1, in order to establish QE, it suffices to prove that V{A) = 0 
for every pseudo-differential operator A of order 0 satisfying E{A) = 0. 

We provide the proof of that fact in the Heisenberg flat case, that is, when —Asr = Rh^h = 
^hRh- For simplicity of notation, we drop the index H in what follows. 

Remark 8. In the Heisenberg flat case, the microlocal factorization of —Asr is global. This fact 
avoids many technicalities. In [6] where the complete proof is done, we have to use microlocal charts 
where the quantum normal form of Theorem 6 is valid. We have also to take care of remainder 
terms, since in the general case we do not have an exact factorization, and handling the additional 
terms raises additional difhculties. 

Let A be a pseudo-differential operator of order 0 satisfying E(A) = 0. Let a = ap{A) be its 
principal symbol. 


Preparation of A. We claim that we can modify A without changing a|s, so as to assume that 
[A, fl] = 0. 

Indeed, it suffices to average with respect to H, as follows. Setting Ag = exp(isSl)Aexp(—isO), 
for s G R, we have, by the Egorov theorem, ap{As) = ao exp(t(Tp(H)). Since (Tp(H) = -|- 

hy)/|hz| = 0 along E, it follows that (Jp{As)\y, = a. Now, setting A = ^ f^^Agds, we have 
‘^p(4)|e = a. Using that exp(2i7rU) = id and that ^Ag = A^], we infer that [U, A] = 0. Then 

we replace A with A. 


Averaging with respect to R. Let us now set At = exp{—itR)Aexp{itR), for t G M.. By the 
Egorov theorem, we have at = ap{At) = a o exp{thz), where hz is the Hamiltonian vector field 
generated by the Hamiltonian function hz- For T > 0, we define Ap = y At dt. To prove that 
U(A) = 0, it suffices to prove that 

U(A-At) = 0, (3) 

and that 


lim V (At) = 0. 

T^+oo 


(4) 


The next two lemmas are devoted to prove (3) and (4). Ergodicity will be used to prove (4). 


Lemma 1. We have U(A — At) = 0, for every t G R. As a consequence, we have U(A — Ap) = 0, 
for every T > 0. 

Proof. Since ^At = i[At,R], it suffices (by Cauchy-Schwarz) to prove that V{[At,R]) = 0. Using 
the Jacobi identity, and the fact that [i?, U] = 0, we compute 


j^[At,n] = z[[At,R]M] 


i[R,[n,At]] + i[At,[R, U]] 


-iad(i?).[At,U], 


and since [A,U] = 0 by construction, we infer that [At,U] = e U] = 0. 

Now, for every integer n, we have 


(: ^n) — 4^n) 4^n) 

— T H” T ; ^n) beCaUSe (pfi — ~ /S-sR^n 

'^n '^n '^n 

= 0 
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because [i?, O] = 0 and = 0. Hence V{A — At) = 0 for every t. 

Using the Fubini theorem and the Jensen inequality, we easily infer that V{A — At) =0. □ 

Lemma 2. We have lim Vi At) = 0. 

T^+oo 

Proof. We have U(Ht) ^ with apiATA^) = |arp = |y atdt\^. By the microlocal 

Weyl law (Theorem 4), we have 

E{AtA*t) = ( \aTiq, ±ag{q))\^ dv, 

^ Jm 


but, using the Reeb flow T^t, 

1 1 

f^iq) = aTiq,±agiq)) = - J atiq,±agiq)) dt = - J a(7^t(q), ±ag(7^^(g))) dt 

By ergodicity of the Reeb flow, and using the von Neumann ergodic theorem, we have 

/t ±<^giq)) dv = 2EiA) = 0 , 

T^+oo 

and the result follows, under the additional assumption that a is even with respect to S. If the 
eigenfunctions are real-valued then this assumption can actually be dropped (see [6]). □ 

4 Perspectives and open problems 

Reeb flow and spectral theory. In addition to QE, there are several other relationships be¬ 
tween the Reeb flow and the spectral theory of sR Laplacians in the 3D contact case. For example, 
in the case where M = S^, the Popp volume is a spectral invariant, which coincides with the 
inverse of the Arnold asymptotic linking number of the Reeb orbits. Moreover, using a global 
normal form, one can figure out a close relationship between the spectrum of an elliptic Toeplitz 
operator that is the quantization of the Reeb flow and the spectrum of As_r. These arguments 
lead to conjecture that the periods of the periodic orbits of the Reeb flow are spectral invariants 
of AsR. 

Measures in sR geometry. A general problem consists of showing the existence of the (micro)- 
local Weyl measures (defined in Definition 2), and of identifying it. 

In contrast to Riemannian geometry, in sR geometry there are several possible choices of in¬ 
trinsic volume forms, for instance: the standard Hausdorff volume (which is defined with arbitrary 
coverings), the spherical Hausdorff volume (which is defined with ball coverings only, with balls 
of constant radius), the Popp volume. The latter one, which is always a smooth volume, was 
introduced and defined in the equiregular case in [17]. Note that for left-invariant (Lie group) sub- 
Riemannian manifolds, the Popp measure and the standard and spherical Hausdorff measures are 
left-invariant, and thus are proportional to the left-Haar measure. In other words they coincide up 
to a multiplying scalar. In particular in the 3D contact case, the Popp measure and the standard 
and spherical Hausdorff measures are proportional, but the constant is not equal to 1, and is not 
known. 

All these measures do not coincide in general (if there is no group structure). More precisely, 
it is proved in [1] that the Popp measure and the spherical Hausdorff measure are proportional 
for regular sR manifolds of dimensions 3 and 4, but in dimension larger than or equal to 5, the 
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spherical Hausdorff measure is not smooth in general and therefore differs from the Popp measure 
(which is always smooth by definition for equiregular distributions). 

The notion of (microlocal, or local) Weyl measure is very natural in the spectral context, but 
seems to be new and has never been studied in sR geometry. It is therefore of interest to compare 
it with other notions of measures and to investigate the following kind of question: does the Weyl 
measure coincide with an appropriate Hausdorff measure? What is its Radon-Nikodym derivative 
with respect to the Popp measure? 

Equiregular sR structures. In Theorem 4, we have obtained a microlocal Weyl law in the 
3D contact case. It is likely that the microlocal Weyl formula can be generalized to equiregular 
sub-Riemannian structures. 

Let {M,D,g) be a sub-Riemannian structure, where M is a compact manifold of dimension 
d. Let Asr be a sub-Riemannian Laplacian on M. We assume that Lie(II) =TM (Hormander’s 
assumption, implying hypoellipticity). We say that the horizontal distribution D is equiregular, if 
the sequence of subbundles (Di)igN* defined by Dq = {0}, Di = D, Di^i = Di + [Di,D], is such 
that the dimension of Di^ dot not depend on q G M, for every integer i. Let r be the smallest 
positive integer such that Dr = TM. We set = Dj-, for j = 0,..., r. 

In the equiregular case, it is known (see [15]) that N{X) ^ CstA'^/^^ with Q = * dim(Di/Di_i) 

(Hausdorff dimension). Besides, as already said in Section 3.2, we already know that, if the 

microlocal Weyl measure exists, then supp(ITa) C SY, = SYi. Actually, we conjecture that 

supp(VPa) = SYr-i, 

and that this result on the support is even valid in general, not only in the equiregular case. 

Contact case in dimension 5. Theorem 2 has been established in the 3D contact case. It 
is natural to investigate the general contact case, however, already in the contact 5D case, the 
situation is much more complicated, and we cannot expect to have a Birkhoff normal form at 
the infinite order as in Theorem 5. Indeed, in dimension 5, the characteristic manifold E is of 
codimension 4, and we have two harmonic oscillators, which may have resonances according to the 
choice of the coefficients of the sR metric. We expect however to be able, in that case, to derive a 

Birkhoff normal form at finite orders, and thanks to a quantized version of that normal form, to 

be able to derive QE under generic assumptions. Indeed, resonances should only occur on sets of 
measure zero in the generic case. 

The role of singular curves. We briefly recall the definition of a singular curve. Let (M, D, g) 
be a sR manifold. A curve q{-) : [0,1] —>■ M is said to be horizontal if it is absolutely continuous 
and if q{t) G for almost every t G [0,1]. Let qo G M he arbitrary. Let fl(qo) be the set of all 

horizontal curves q(-) such that q(0) = qg- It is clear that D(go) is a Banach manifold, modeled on 
L^(0,1;]R™) (with m = rank(D)). Now, let qi G M he another arbitrary point. The set Vl{qQ,qi) 
of all q{-) G fl(go) such that g(l) = qi may not be a manifold. Indeed, defining the end-point 
mapping endq^ : fl(go) -)> M by endg(,(g(-)) = g(l), we have fl(qo,qi) = (endgj~^(gi), and the 
mapping endg^ need not be a submersion. By definition, a singular curve q(-) is an horizontal curve 
that is a critical point of endg(Q), or, in other words, such that the differential dendg(o) is not of 
full rank. 

There are many clues leading one to think that singular curves may have a strong impact 
on the asymptotic spectral properties of the sR Laplacians. For instance, we expect that they 
have an impact on the Schwartz kernel in the representation through Feynman integrals. This 
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representation has the form 


exp(zt/iA)/(a;) = f d. 

J g(')Gf2(x) 

and if we fix the terminal points of the paths, then we “disintegrate” the measure to obtain the 
Schwartz kernel. We expect that singular curves cause a singularity in the kernel, which differs 
from the usual one. In order to understand such issues, we provide hereafter several examples. 


Martinet case. We consider a smooth compact connected manifold M of dimension 3, and a 
generic horizontal distribution D of rank 2, defined hy D = kera, with a a nontrivial one-form on 
M such that a A da vanishes transversally on a surface S of dimension 2. This surface is called 
the Martinet surface. Let g be a Riemannian metric on D. 

The local model in near 0 is given hy a = dz — dy. Then D is locally spanned by the 
two vector fields X = dx and Y = dy + x^dz- Note that S' = {a; = 0} and that [X, T] = 2xdz 
vanishes along S. A Lie bracket of length two, [A, [A, A]] = 29^, is required to generate the missing 
direction, along S. Outside of S, the distribution D is of contact type. 

Note that the characteristic manifold S is not symplectic. This is equivalent to the existence 
of nontrivial singular curves. The Martinet case is indeed a well known sR model in which there 
are nontrivial singular curves, that are minimizing, and locally foliate the Martinet surface S (see 

[17]). 

We assume that L = DOTS is a line bundle over S (otherwise, in the general case, L may have 
generic singularities): hence S is either a torus or a Klein bottle. Let dP be the Popp volume on 
M \ S. Near S, we have dP ~ dv 0 |^|, where v is a smooth measure on S and </> = 0 is a local 
equation of S. 

Given any smooth measure /i on M, we consider the corresponding sR Laplacian, and some 
eigenbasis of it. We are able to obtain the following local Weyl law (see [7]). 


Proposition 1. For every 


continuous function f on M, we have 


X! / X^logX, 


An<A 


/M 


32 


as \ ^ -boo. In particular, N{\) ^ log A. 

Corollary 1. There exists a density-one subsequence (nj)jgN* sueh that 


lim f f\(j)n fdy = 0, 

for every continuous function f on M, vanishing near S. 

Note that Proposition 1 cannot be obtained from the asymptotic estimates of the heat kernel 
along the diagonal, because those estimates are not uniform: along the diagonal, the asymptotics 
is in on S, and in C{q)t~^ outside of S, with C{q) bounded but not integrable near S. The 
proof of Proposition 1 consists of first computing the asymptotic behavior of Tr(/e“‘^) as t —>■ 0+ 
in the flat Martinet case, by rescaling. Then, using the locality of the heat kernel expansion along 
the diagonal, the result is obtained in the locally flat case. Using Dirichlet-Dirichlet bracketing, 
the result is then established in small cubes with a flat metric, and is finally extended using local 
quasi-isometries and minimax considerations. 

It might be expected that the singular flow plays an important role in the asymptotic spectral 
properties of the sR Laplacian. In particular, it is tempting to think that a good assumption is 
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the ergodicity of the abnormal geodesics (lifts of the singular curves in the cotangent space) in 
the characteristic distribution kerwis, at least, under the assumption that there indeed exists a 
characteristic vector field leaving v invariant. 

Relationship between singular curves and magnetic fields. This paragraph is inspired by 
the paper [16]. Let us consider a manifold N of dimension 2, endowed with a flat Riemannian 
metric and with its canonical Riemannian volume. We consider the Laplacian 

A = -{d^ - ad^f - {dy - bd^f, 

with z G 'R.I2'k'L and A = adx + bdy the magnetic potential on N. Setting X = dx — adz and 
Y = dy — bdz, we have [X, Y] = B dz where dA = B dx A dy is the magnetic field. We assume that 
B vanishes on a closed curve L with a non-zero differential. Then a = a dx + b dy + dz defines 
a Martinet distribution on M = N x R/2 ttZ, and the Martinet surface is S' = L x R/27rZ. We 
have, then, Og = ^{adx + bdy + dz) and ix = -pj^lds dzj, where s is the arc-length along T. The 
previous expressions are still valid if the metric is not fiat. The characteristic direction is given 
by d.( 7 (s)) ds -b dz = 0 if r is parametrized by the arc-length. It is always possible to choose a 
gauge A near T such that ^.( 7 ) = Aq is constant. Using this gauge, the measure zz is invariant 
under the vector field ||dR||(ds — Aodz). The corresponding dynamics is ergodic if 27rzlo/length(r) 
is irrational. 

Another case of interest is the quasi-contact case in dimension 4. It is related to magnetic fields 
in dimension 3. In this case, there exist some nontrivial singular curves that correspond to lines 
of a magnetic field. It seems that the singularities can only occur at isolated points. We might 
then expect to have QE if there is only one such point, and otherwise a tunnel effect might occur 
between two different points. 

Almost-Riemannian geometry. For the sake of simplicity, we restrict ourselves to almost- 
Riemannian (aR) structures on surfaces. Let M be a smooth compact connected manifold of 
dimension 2. An aR structure on M is locally given by two vector fields X and Y that generate, 
outside of a singular curve 5, the tangent space of the surface M. We assume that, on S, the 
vector fields X, Y and [AT, T] generate TM. The simplest example is the so-called Grushin case 
and is given locally hy X = dx and Y = xdy. Outside of S', the vector fields X and Y define the 
Riemannian metric g = + Vy, where [vx^vy) is dual to (A,T). The volume y,aR associated 

with g is singular along S. In the Grushin case, we have g = dx^ + {dy/x)^ and dy,aR = \dx dy\/\x\. 

Let dfi be any volume form on M. We consider the aR Laplacian AaR = —X*X — Y*Y, where 
the adjoints are taken with respect to /r. It follows from Remark 7 that the aR Laplacians associated 
with two different volumes are unitary equivalent up to a potential of the order of d{-,S)~^. If 
M = fJ-aR then the resulting operator is singular on S, however it is essentially selfadjoint on M\S 
(see [3]). 

In this Grushin case, we prove in [7] a result that is similar to the Martinet case, namely, that 
for every continuous function / on M, we have 

as A —>• -boo (with, as well, Corollary I as a consequence). 

Remark 9. In any of the above singular cases (Martinet, almost-Riemannian), one may think of 
applying the desingularization procedure of [19]. For instance, the Martinet case can be seen as the 
projection of the so-called Engel case, which refers to an equiregular horizontal distribution of rank 
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2 in dimension 4; the Grushin case can be seen as the projection of the 3D contact case. However, 
in view of obtaining a microlocal Weyl law, it is not clear whether or not this desingularization 
can be used in a relevant way, because when projecting one has to use integrals of the Schwartz 
kernel outside of the diagonal, that are not known. 

Remark 10. In order to study sR structures in which the rank of the horizontal distribution D 
may not be constant, a more general definition of a sR structure has been introduced (although 
not exactly with those terms) in [2] . 

A relative tangent bundle on a manifold M is a couple where tt : R —>■ M is a smooth 

fibration and ^ H TM is a smooth bundle morphism over M. The horizontal distribution is 
defined hy D = ^{H) C TM. A vector field X on M is said to be horizontal if there exists a section 
u oi H such that X = ^{u). An horizontal curve q{-) : [0, 1] —>■ M is an absolutely continuous curve 
for which there exists a section u of H such that q{t) = for almost every t G [0,1]. 

A sR structure (with possibly varying rank) is then more generally defined by a quadruple 
(M, H, g), where M is a manifold, (H, is a relative tangent bundle on M, and g is a Riemannian 
metric on H. This definition may be used in order to provide a unifying context for usual sR and 
aR geometries. 

Microlocal properties of SR wave equations. It is interesting to investigate properties of 
wave equations associated with a sR Laplacian. One of the first questions that are arising is the 
question of microlocalization. More precisely, let us assume that the wave front of Cauchy data 
(uqjUi) is a subset of a cone K with compact base contained in an open cone C. Does there 
exist e(K, C) > 0 (not depending on the Cauchy data) such that the wave front of the solution 
is contained in C for |t| ^ e? This question is a priori far from trivial because the sR cut-locus 
and the sR conjugate locus of a point always contain this point in their closure. Another natural 
(simpler) question is to establish finite propagation of waves in the sR context. 

Quantum limits. The decomposition and the statement on /Iq given in Theorem 3 are actually 
valid for any sR Laplacian, not only in the 3D contact case. A stimulating open question is to 
establish invariance properties of the measure (3oo for more general sR geometries. 
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